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Abstract 

In this work we consider the Navier-Stokes problem modified by the 
absorption term |u|''^^u, where cr > 1, which is introduced in the mo- 
mentum equation. For this new problem, we prove the existence of weak 
solutions for any dimension N > 2 and its uniqueness for A'' = 2. Then 
we prove that, for zero body forces, the weak solutions extinct in a finite 
time if 1 < (J < 2, exponentially decay in time if cr = 2 and decay with a 
power-time rate if ct > 2. We prove also that for a general non-zero body 
forces, the weak solutions exponentially decay in time for any cr > 1. In 
the special case of a suitable forces field which vanishes at some instant, 
we prove that the weak solutions extinct at the same instant provided 
1< cr < 2. 

MSG: 35Q30, 76D03, 35B40. 

Keywords: modified Navier-Stokes, absorption, existence, uniqueness, extinc- 
tion in time, power decay, exponential decay. 

1 Introduction 



1.1 Basic equations 

From the basic principles of Fluid Mechanics, it is well known that, in isothermal 
motions of incompressible fluids, the velocity field and pressure are determined 
from: 

• the incompressibility condition 

divu = 0; (1.1) 
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• the conservation of mass 

|^+div(pu) = 0; (1.2) 

• the conservation of momentum 

p(j^ + (u- V)u^ =pf+ divS. (1.3) 

The notation used in l|l.ip - (|1.3p is well known: u is the velocity field, p is the 
pressure, p is the density and f is the forcing term. For Newtonian fiuids, the 
stress tensor S is given by the Stokes law 

S = -pI + 2^D, D = i(Vu + Vu^), (1.4) 

where ^ is the dynamical (constant) viscosity and D is the strain tensor. For ho- 
mogeneous fiuids, the density is regarded as constant. Therefore, we can replace 
the continuity equation (|1.2p by its incompressible form l|l.ip . In consequence, 
when these simplifying features are present, the equations for a linearly viscous, 
homogeneous, incompressible fiuid reduce to the following system 

+ (u-V)u = f- -Vp + z^Au (1.5) 

divu = 0. (1.6) 

Here ly = fi/p is a. constant corresponding to the denominated kinematics vis- 
cosity and, because the fiuid is homogeneous, p is a constant corresponding 
to the initial density. In the literature, the system of equations (|1.5p - (|1.6p is 
known as the incompressible Navier-Stokes equations, or only the Navier-Stokes 
equations (NSE). System l|1.5p - l|1.6p must be supplemented with boundary con- 
ditions characterizing the fiow on the boundary of the domain occupied by the 
fiuid and by initial conditions determining the initial state of the fiow at the 
beginning of the time interval. The question of initial conditions is immediately 
understood from the physical point of view, but as for the boundary conditions 
is much more delicate and would require a detailed discussion. See Serrin [21] 
for a detailed derivation of the Navier-Stokes equations and for the explanation 
of the initial and possible different boundary conditions. 

1.2 The modified Navier-Stokes problem 

Let us consider a general cylinder 

Ot := X (0, T) C X E+ , with Ft : dn x (0, T) , 

where is a bounded domain with a compact boundary dQ. Although the 
dimensions of physical interest are N = 2 and = 3, here we shall consider a 
general dimension N > 2. In this work we consider the following modified NSE 

-^ + (u-V)u = f- Vp + z/Au-ajur^^u in Qt , (1.7) 
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divu = in Qt 



(1.8) 



supplemented by the initial and boundary conditions 



u = uo when t — Q in 



(1.9) 



u = on Tt ■ 



(1.10) 



In l|1.7p . a and a are positive constants with a > 1. Notice that, for the sake of 
simplicity, we have assumed p = 1 in l|1.7|) . In the sequel (|1.7p - l|1.10p shall be 
denominated as the modified NS problem. The motivation for the consideration 
of the absorption term |u|'^~^u in l|1.7p is purely mathematical and goes back to 
the works of Benilan et al. [7], Diaz and Herrero [TT], and Bernis 0|9]. There, 
was studied the importance of a similar absorption term to prove qualitative 
properties related with compact supported solutions, or solutions which exhibit 
finite speed of propagations, or which extinct in time, for different initial and or 
boundary value problems. A possible physical justification for the absorption 
term in (|1.7p . is the consideration, in the momentum equation, of a forces field 



where f is a given forces field. Notice that such forces field depends, in a sub- 
linear way, on the own velocity u. In a certain sense, such forces field may 
be considered, from the physical point of view, as a feedback field. Therefore, 
l|1.7p can be considered as the equation of motion for a certain fiuid. In this 
way, the purpose of this work is to study the response of such a fiuid motion 
undergone to a forces field satisfying (|l.lip . Specifically, we want to know if 
such a forces field is responsible for stopping the fiuid, driven by (|1.7p - (|1.10p . 
This issue addresses us for the important question about the decays of the solu- 
tions of the Navier-Stokes equations. With this respect, it should be remarked 
that these questions have been studied by many authors, among many others, 
by Horgan and Wheeler [15], Wiegner |25], Ames and Payne [T], Schonbek [22] . 
Borchers and Miyakawa [10], Kozono and Ogawa [16], Takahashi [23], Enomoto 
and Shibata [12] and Bae and Jin [5] . In the references [TS^ [1] are established 
the exponential spacial decay for the solutions of the stationary Navier-Stokes 
equations. The other references [Ml [22l SOl SSI O [l2l E] are concerned with 
space and time power decays in different norms for the solutions and its deriva- 
tives of the Navier-Stokes equations supplemented with suitable conditions on 
the initial data or in the forces field. In spite of many work in this field, so far, 
and to the best of our knowledge, there are no results establishing the extinction 
of the solutions of the Navier-Stokes equations in a finite time or in space. 

1.3 Mathematical framework 

Notation. The notation used throughout this text is largely standard in Math- 
ematical Analysis and in particular in Mathematical Fluid Mechanics - see, e.g., 
Lions [18] , Temam [24] or Galdi [13} fT4] . We distinguish vectors from scalars 
by using boldface letters. For functions and function spaces we will use this 



h(u) = i-a\M\''-^M, 



(1.11) 
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distinction as well. The symbol G will denote a generic constant - almost the 
time a positive constant, whose value will not be specified; it can change from 
one inequality to another. The dependence of C on other constants or parame- 
ters will always be clear from the exposition. Sometimes we will use subscripted 
letters attached to C to relate a constant with the result it is derived from. In 
this article, the notation Vt stands always for a domain, i.e., a connected open 
subset of M^, whose compact boundary is denoted by dO,. 

Function spaces. Let 1 < p < oo and 51 C be a domain. We shall use the 
classical Lebesgue spaces \J'(yi), whose norm is denoted by || • ||lp(o)- For any 
nonnegative fc, W'^'P(ll) denotes the Sobolev space of all functions u S l-Pip) 
such that the weak derivatives D"u exist, in the generaHzed sense, and are in 
L''(r2) for any multi-index a such that Q <\a\ < k. The norm in W'^'P(ri) is de- 
noted by II • |lw''.p(n)- The associated trace spaces are denoted by W'=-i/P'P(9n). 
Given T > and a Banach space X, LP(0,r;X) and W^'P{Q,T;X) denote the 
usual Lebesgue and Sobolev spaces used in evolutive problems, with norms de- 
noted by II • ||lp(o,T;X) and || • ||w'' p(o.T:X)- The corresponding spaces of vector- 
valued functions are denoted by boldface letters. All these spaces are Banach 
spaces and the Hilbert framework corresponds to p = 2. In the last case, we use 
the abbreviations W^-^ ^ jjfe ^'=-1/2,2 ^ h'=-i/2. 

Auxiliary results. Throughout this text we will make reference, more than 
once, to the following inequalities: 

(1) Algebraic inequality - for every a, /3 G M and every A, B >Q, 

A^'B^ <{A + By''+^- (1.12) 

(2) Young's inequality - for every a, 6 > 0, e > and 1 < p, q < 00 such that 
l/p+l/q=l, 

ab<eaP + C{e)b\ C{e) = {epy^^/Pq-^ . 
li p = q = 2, this is known as Cauchy's inequality. 

(3) Holder's inequality - for every u £ L^'(fJ), v e L''{n), with 1 < p, q < co 
such that l/p + l/q = 1, 

uvdyi < \\u\\i^p^i2)\\v\\Li{n)- 

An important result that will be used in the sequel is the famous Gagliardo- 
Nirenberg-Sobolev inequality. 

Lemma 1.1 Let be a domain ofM.^, N > 1, with a compact boundary d^l. 
Assume that u € Wj'''(ri). Then, for every fixed number r > 1 there exists a 
constant C depending only on N, p, r such that 

hllL,(o) <C^I|Vii||^.(j,)h||J,7fo). (1-13) 
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where p, q > 1, are linked by 



1 1 



111 



r q J \N p r 



(1.14) 



and their admissible range is: 
1. //7V=1, 9G[r,oo], 0e 



0, 



p 



2. Ifp < N,qe 



Np 



p+r{p—l) 

Np 
N~p 



ifr >j^andqe 



Np 
N~p 



«/r<^,0G[O,l] 



C=[iN-l)p/iN-p)]' 



3. If p > N > 1, q e [r,oo), 9 e 
l)/N,l + {p-l)pry. 



Np 



Np+r{p-N) 



and C = max{q{N 



When 9 — 1, (|1.13|) is known as the Sobolev inequality and, in this case, if 
q = p ^ 2, then l|1.13p is usuaUy denominated as the Poincare inequality. This 
result is valid whether the domain fl is bounded or not and notice the constant 
C does not depend on fl. See the proof in Ladyzhenskaya et al. [iTl p. 62]. The 
following generalization of l|1.13p to higher-order derivatives is also possible (see 
Nirenberg [19]). 



^.|L,(n)<^l|D'-llL.(0)l 



Il'-(o) 



(1.15) 



where additionally < j < fc and A: > 1, and now 



N 



k 1 

N^p 



The results written above can be easily generalized for vector-valued func- 
tions. Moreover, for time-dependent functions u{x,t), they still hold for a. a. 

[o,T]. 

Other very important auxihary result which will be of the utmost importance 
to handle the absorption term is written in the following lemma. 

Lemma 1.2 For all p G (l,oo) and S > 0, there exist constants Ci and C2, 
depending on p and N, such that for all ^, 77 G E^, iV > 1, 



\\^r'^~\vr'v\<ci\^~v\'-'m 



\v\ 



\p-2+S 



and 



(ler-'^ - ivr'v) • (e - ^) > c^2ie - m + \v\ 



.p-2-S 



(1.16) 



(1.17) 



See Barret and Liu [6] for the proof and also the references cited therein for 
other forms of ifTTTe)) and (fTTZl) . 
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2 Weak formulation 



In this section, we establish existence and uniqueness results for the modified 
NS problem pT7l) - (|1.10p . Notice that, in the limit case a = in (fLTj) . we fall 
in the classical NS system (|1.5p - l|1.6p . In this case, it is well known that the 
corresponding problem has a weak solution which is unique only if iV = 2 (see 
e.g. Lions [18] and Galdi [H]). In order to define the notion of a weak solution 
to the problem l|1.7p - (|1.10p . let us introduce the function spaces largely used in 
Mathematical Fluid Mechanics: 

V:={veCS°(f^) :divv = 0}; (2.18) 

H := closure of V in L^(f7) ; (2.19) 

V, := closure of V in H" (17), s>l; (2.20) 

where, for simplicity, we can assume s as the smaller integer not lesser than 
N/2, to avoid the compHcated Sobolev spaces with s non- integer. H is endowed 
with the L^(r2) inner product and norm, and, for any positive integer s, is 
endowed with the inner product 




and with the associated norm. From (|1.13p . we see that this norm is equivalent 
to the HQ(ri) norm and, in consequence, we have the (isometrically isomorphic) 
identification = H~'''(f7). H and are Hilbert spaces and, for s = 1, which 
happens when the dimension is = 2, we simply denote, as usual, Vi by V. 
In this case, the V inner product reads in the usual form 

(u, v)v= / VuiVvdx. 
Jn 

From the theory of distributions, we know that 

V H = H' V' V;, , s > 1 . (2.21) 

Moreover the compact imbedding Hj(r2) ^ L^(17) implies that the imbedding 
Vs ^ H is also compact for any s > 1 (see Lions [TSt pp. 67-77]). For the theory 
of these function spaces, see Lions [18], Temam |24| and Galdi [13]. The notion 
of weak solution for the problem (|1.7p - (|l.lQp follows in a standard manner. 

Definition 2.1 Let Q be a bounded domain in ii^, N >2. A vector field u is 
a weak solution of the problem ^.7]) - lll.lO\) . if a > 1 and: 

1. ue L^{0,T;V}r\L''{QT)nL°°{0,T;H); 

2. tt(-,0) ~ uq a.e. in fl; 
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3. For every ve VD (Q) n X''(f7) and for a.a. t e (0,r), 




(2.22) 



+ a \u{t)\''^'^u{t)- vdx= / f{t)-vdx. 




Condition u G L^(0,r;V) expresses, in a certain sense, the incompressibil- 
ity condition (|1.8p and the boundary condition l|1.10p . whereas condition u e 
L°°(0,T; H), a priori, does not seem to be strictly necessary or else to restrict 
the class of admissible weak solutions. Condition u G L"{Qt) is a natural re- 
quirement to deal with the absorption term |u|'^^^u. Condition u(-,0) = uq 
a.e. in fi, should be interpreted in the sense that u is L^(r2) weakly continuous 
for t = 0, i.e. 



For and cr < 4, l|2.22p holds for every v G V, because, due to Sobolev's 
inequality, H\n) ^ V\n) for p < 4 and, in consequence, V n L^(ri) = V. 
Definition 12.11 is silent about the initial data uq and the forces field f. But, this 
will be clear when we bellow establish the existence result. 

2.1 Existence result 

To prove the existence of a weak solution of the modified NS problem l|1.7p - 
l|1.10P , we will adapt the same arguments used to prove existence for the classical 
NS problem (see, e.g. Lions [18], Temam [51], Galdi [H]). However, it is worth 
to notice that, in (|1.7p . additionally to the usual nonlinear term for the classical 
NS equations, (u-V)u, we have another one, the absorption term a|u|°'~^u. We 
shall adapt the proof for the classical NS problem in any dimension N >2 given 
in Lions jlHl pp. 75-77]. 

Theorem 2.1 Assume that f € L^(0,T; F') and Uo & H. Then, there exists, 
at least, a weak solution of the modified NS problem ii. 7)) -( f7TiQ)) in the sense of 
Definition \2.1\ 

PROOF. 1. Existence of approximate solutions. Let s > 1 be the smaller integer 
not lesser than A^/2. We may assume that s > 1, because for s = 1, V,, = V 
and the proof would follow in the same manner, but even simpler. We consider 
the basis {vfc}j,g]y[ of Vs, given by the (non-zero) solutions \j of the following 
spectral problem associated to the eigenvalues \j > 0: 



Since Vs ^ V ^ H, {vfc}j.gjsj can be chosen as being an orthonormal basis 
in H. Let us consider also the corresponding m-dimensional space, say V™, 






D"vj • D"v3 dx = Aj / Vj • dx V (/? G Vs . 
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d 

dt 



spanned by vi, . . . , v„i. For each to G N, we search an approximate solution 
Um of l|2.22p in the form 

m 

Um = ^ Ckm{t)vk , (2.23) 
k=l 

where Vfe S V™ and Ckmit) are the functions we look for. These functions are 
found by solving the following system of ordinary differential equations obtained 
from l(232l) : 

/ u„i(i) -Vkdx + iy / V(u„i(i)) : V Vfcdx + / (u„(t)-V)u„j(t) • Vfedx 

f f (2.24) 

+a |um(i)|'^ u„(t)-Vfc(ix= / f(t)-v/c(ix; 

Cfcm(O) = / uo„ • Vfcdx; (2.25) 

for k ~ 1, . . . ,m and where uom = Um(0) £ V™ is such that 

uom — > uo strongly in H as to — > oo. (2.26) 

From the elementary theory of ordinary differential equations, problem l|2.23p - 
l|2.25p has a unique solution € C^([0,rm]), for some small interval of time 

[o,r„] c [o,T]. 

2. A priori estimates I. We multiply l|2.24p by Ckm{t) and add these equations 
from k = 1 to k = m. Then, according to Lions [TSt Lemme 1-6.5] (see also 
Temam [24l Lemma II-1.3]), we obtain 

^;^l|u™(t)|!^(f,) +i^||Vu™(t)||2,(^j^ +a||u™(i)||£„(j,) 
^ (2.27) 
u,„(t) • f{t) dx . 

Integrating (|2.27p from to i < T^, using Schwarz's and Cauchy's (with a 
suitable e > 0) inequalities, we achieve to 

l|umWllL2(j2) +^ / l|Vu™(s)||^2(j-jjds + 2a / ||u™(s)||^„(j^)ds < 

^ (2.28) 



for t < Trn. Since ||uom IlL=^(n) < ||uo||L2(n), the assumptions uq G H and 
f G L'^(0, T; V') justify that the left-hand side of l|2.28p is finite. In particular, 
it follows that |cfcm(t)|^ < oo for all fc = 1, . . . , to. In consequence, by standard 
results on ordinary differential equations, we get that T„j = T for all to G N, 
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otherwise |cfem(i)| ^ oo as t — > T„i. Moreover, from (|2.28p and once that 
a > 0, we obtain 

sup ||u„,(i)||^(o) < I|uo|Il2(o) + -llfll^(oT-v') ' (2-29) 

which impUes that 

Um remains bounded in L°°(0,r;H). (2.30) 
On the other hand, if we replace i by T in l|2.28p . we obtain 



< 



(2.31) 



I"0|Il2(O) + -||f|lL2(0,T;V') 



1 
V 

This estimate enables us to say that 

Um remains bounded in L^(0,r; V) and in Ij"{Qt). (2.32) 

Moreover, 

|um|'^~^Um remains bounded in {Qt)- (2.33) 

3. A priori estimates II. Let us consider the orthogonal projection P,„ : H 
V", Pm{-a) = X;r=i L " • '^fe dyi^fk. We thus obtain from l(2^ 

= (l/Au„) - Pra ((U„ • V)U„0 - Pra (|u„|"-2u„) + P„, (f) , (2.34) 

where we have used Pm(u„j) = u„i by virtue of (|2.23p . Using l|2.29p and (|2.3ip . 
and the special choice of the basis of V™, we deduce, arguing as in Lions [l8t 
1-6.4.3], that the sequences Pm {vl^vLm) and Pm ((um ■ V)um) are bounded in 
L^(0,T;Vg). Moreover, from hypothesis, one clearly has f S L^(0, T;Vg). Then, 
from l(2:34| . 

remains bounded in L2(0, T; V'J + 1,"' {Qt). (2.35) 

4. Passing to the limit. From l(230| . l(233| and (|235l) . there exist 

functions u and v, and there exists a subsequence, which we still denote by u^, 
such that 

u,„ ^ u weak-star in L°°(0,T;H) as ?7i — > 00 , (2.36) 

Um ^ u weakly in L^(0,r;V) as m cxo , (2.37) 
u,„ — » u weakly in Ij" {Qt) as m ^ 00 , (2.38) 
|um|'^~^u„i — > V weakly in \f {Qt) as m ^ 00 (2.39) 



9 



and 

~Wr'^^. weakly in h'^{Q,T]\[^) + 1."' {Qt) as m ^ oo . (2.40) 

Then, due to (|2.2ip . (|2.37p and l|2.40p . and according to a well known compact- 
ness result (see Lions [TSl Theoreme 1-5.1]), 

Um — > u strongly in L^(0,T;H) as m ^ oo . (2-41) 

Now, we multiply (j2.24p by V' G C^([0, T]), with '0(T) = 0, and then we integrate 
the resulting equations from to T. We thus obtain 

- / / Um(t) • Vfe i^'it) dxdt + v ( [ Vu,„(i) : V V(0 dxdt+ 
Jq Jn Jo Jn 

{urn{t)-V)u^{t)-Yk^p{t)dxdt + a f f \u^{t)\''-^u^{t)-xk^p{t)dxdt 

jQ Jn 

= 11 m ■ Vfe ij{t) dxdt + [ u,nO ■ Vfe V(0) dx . 
Jq Jn Jn 

(2.42) 

Now, we can pass to the limit in l|2.42p by using l|2.37p in the Hnear terms (|2.42P i 
and 1(2^ 9. l|237l) along with ((2^ in the nonlinear term ([231)3 and l[236l) 
in the term (|2.42P k (see Lions [HI 1-6.4]). With respect to the term l|2.42p .d. we 
notice that by considering a new subsequence, we may assume that u„ u a.e. 
in Qt- This and (|2.39p prove that v = |u|°'~^u, and we can pass to the Hmit in 
this term as well. Finally, using linear and continuity arguments, we can prove 
that u satisfies l|2.22p in the distribution sense of CJ^(0,T). Moreover, using 
standard arguments, we can prove that u satisfies 2. of Definition 12.11 See e.g. 
Temam [H pp. 288-289] for the details. □ 

Remarks 2.1 1. This existence result still holds for unbounded domains fl as 
far as the used Sobolev type inequalities hold. The main difference lies in the 
fact that the imbedding V ^ H is no longer compact. But this difficulty is 
overcame in a standard manner (see e.g. Temam 124\ Remark IIT3.2]). 
2. From what we have said at the beginning of this section, in particular I12.21\) . 
the assumption f <E L^(0,T; V') can be replaced by f £ L^{0,T; ir^{fl)) . On 
the other hand, with some minor modifications in the proof, we can assume that 
/G L-^(0, T; Z^(ri)) for some some p > 1. 

2.2 Energy relation 

In this section we shall establish the energy relations satisfied by the weak 
solutions of the modified NS problem l|1.7p - l|1.10p . Before we estabHsh the main 
result, let us define the (kinetic) energy associated with the problem (|1.7p - (|1.10p : 

m--=l\Ht)\\l.^n)- (2-43) 
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As one should expect, at least from the physical point of view, any solution 
of problem (|1.7p - l|1.10|) should be such that the associated kinetic energy at a 
certain time t (considering f = 0) is equal to E{s), for some < s < t, minus 
the energy dissipated by viscosity in [s,t], 

V ||Vu(T)||^2(f^)dr (enstrophy), 

and minus the energy dissipated by a sink associated to a in [s,t], 

"/ llu(T)||L.(f^)dT (absorption). 

J s 

However, even for the classical NS problem, weak solutions obey to such an 
equality only if — 2. For iV > 3, it is only possible to prove they satisfy to 
an energy inequality (see e.g. Galdi [Mj). On the other hand, from the above 
considerations, we should expect 

uo G H ^ E{Q) <oo^ E{t) < E{0) < oo for alH > . (2.44) 

In the following result we shall see the weak solutions of the modified NS problem 
l|1.7p - l|1.10p satisfy to an energy inequality, regardless the domain dimension. 

Theorem 2.2 Assume that f G L^(0,T; V) and Uq £ H. Let u he a weak 
solution of the problem Ill.7\} - lll.l0]) in the sense of Definition \2.1[ Then this 
solution satisfies to 

^4ll«WII^(o)+'^l|Vu(t)||2,(j^) +a||M(t)||^„(o) < 

r (2-45) 
/ u{t)-f{t)dx 

for a. a. t G [0,r], and 




(2.46) 



for a.a. te [0,T]. 

PROOF. To prove ([TISl) we take the Hmit inf, as m ^ oo, of ^Mil- Then from 
1^261) . ([236]) - ((2381) and a classical property of weak limits (see Galdi [I^), the 
result follows. For (|2.46p . we start by integrating (|2.27p between and t < T. 
We thus obtain 




(2.47) 
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Now, we take the limit inf, as m -> oo, of (|2:47|) . Then from ([236)) - l(238)) 
and the same classical property of weak limits, it follows l|2.46p . □ 

Such as for the classical NS problem, we can prove the weak solutions of the 
modified NS problem (fLTil - lfrTOl) belonging to L^{0,T]L'^{n)) satisfy to the 
energy equality 



for a. a. t e [0, T]. The proof can be adapted from Galdi |14[ Theorem 4.1]. The 
only difference lies in the absorption term which can be handle such as in the 
proof of Theorem 12.11 Moreover, if A'' = 2, every weak solution of the modified 
NS problem satisfies to the energy equality (|2.48|) . Indeed for iV = 2 any such 
weak solution u satisfies to the so-called Ladyzhenskaya inequality 



This is no longer valid for N > 3 and therefore, such as for the classical Navier- 
Stokes problem, the question of whether a weak solution of the modified NS 
problem (|1.7p - (|1.10p obeys the energy equality l|2.48p remains open. 

2.3 Uniqueness result 

In the mathematical theory of the NS equations another important open prob- 
lem is the uniqueness of weak solutions for > 3. It is possible to show 
uniqueness for all time under additional conditions on the admissible weak solu- 
tions, or in a small interval of time by assuming small data compared to viscosity 
(see e.g. Galdi [H]). However, for N = 2 the mathematical theory is complete 
and the existence result is now classical. In consequence, we will show the 2-D 
modified NS problem (|1.7p - (|1.10p inherits this property. The proof is again an 
adaptation of the classical result (see e.g. Lions [H]) and the crucial part lies 
in the absorption term. 

Theorem 2.3 Assume that fe L2(0,r; V') and e H. Let V and w be two 
weak solutions of the modified NS problem fi.7)) - f7775)) in the sense of Defini- 
tion \2.1[ Then v a.e. in Qt- 

PROOF. Arguing as in Temam [211 Theorem III-3.2], we get from (|2?22|) the 
following relation for u = v — w: 





(2.48) 





+ 2t.|lVu(t)||2, +/i =/2, 



(2.49) 
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where 



12-.= 2 [ [(w(t)-V)w(t) - (v(t)-V)v(t)] •u(t)dx. 



In the term Ii, we use (|1.17p with ^ = v, = w, p = cr and (5 = 0, to prove that 
h > 0. On the other hand, it can be proved, in a standard manner, that 



for some positive constant. Then, using Poincare's inequahty, we obtain, from 
1(2^ and ((2?50l) . the following relation 



Integrating (|2.5ip . using l|2.3ip for w, and known that u(0) = 0, we prove that 
V = w a.e. in Qt-O 

For TV > 3 and analogously as to the classical NS problem, one can proves the 
weak solution of the modified NS problem, satisfying l|2.46p . is unique in the 
class L2(0, T; V) n L°°(0, T; H) n L^^Qt) n L"(0, T; L''(rj)) with 2/s + N/r<l 
and for r > (see Lions |lHl Theoreme 1-6.9]). In particular the uniqueness 
holds in the class L^iO, T; V) n L°°(0, T; H) n L'^(Qt) for a>N + 2. 

3 Asymptotic stability 

In this section we shall study the behavior in time of the weak solutions of the 
modified NS problem l|1.7p - l|1.10p . Here, we shall assume that T is sufficiently 
large or even let T = cxd. It is worth to recall that, to the best of our knowl- 
edge, the late studies on the asymptotic behavior of the weak solutions of the 
classical NS problem provide only power-time decays (see the references cited in 
Section [L2|) . The first result we present here, estabHshes a temporal qualitative 
property satisfied by the weak solutions of the modified NS problem (|1.7p - (|1.10p . 
which is usually referred to as the extinction in (a finite) time property. 

Theorem 3.1 (Extinction in time) Assume 1 < a < 2 and uq E H, and let 

u be a weak solution of the modified NS problem ^^-iLj^ in the sense of 
Definition \2.1l 

1. If f = a.e. in Qt, then there exists t* > such that u = a.e. in Vl 
and for almost all t > t* . 



\h\ < ^||u(t)||^i(o) + 7l|uWII^(o)l|w(t)||^i(a) , 



(2.50) 



;|l|u(^)ll^(o)<^l|u(^)II^(o)l|w(^)llL^(o)• 



(2.51) 



2. Iff^ 0, but 




(3.52) 
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where tf is a fixed positive time and fi is given by fS. 58\} . then there exists 
a constant eq > such that u ~ a.e. in and for almost all t > tf, 
< e < eo. 



PROOF. First Step. If f a.e. in Qt, we obtain from f^A5\i 

^E{t) + CE2^a (t) < for a.a. t € [0, T] , (3.53) 
where E{t) is the (kinetic) energy defined in (|2.43p . C = min(i^, a) and 

E^.At) l|Vu(0||^(o) + l|u(i)||^(o) • (3-54) 

Using a vectorial version of the GagHardo-Nirenberg-Sobolev inequaUty l|1.13p 
with p = q = 2 and r — a, 

l|u(t)||L^(0) <CGivs||Vu(t)|i^.(f,)||u(0||[^fo) fora.a. <G [0,r], (3.55) 

where Cqns — C{N,a) is the constant resulting from applying i|1.13p and, 
according to (|1.14p . 

= 1- — — ■ (3.56) 

Then, using the algebraic inequality (|1.12p . we obtain 

\\^m\h(n) < Clt,sE2Atr for a.a. t G [0,T] , (3.57) 
where, from p.55p and l|3.56p . 

2(2 -a) 

^-^^+ (2^;)iV + 2. - ^'-''^ 

Then (|3.53p and l|3.57p lead us to the homogeneous ordinary differential inequal- 
ity 

j^E{t) + CE{tY'^' < for a.a. t G [0, T] , (3.59) 
where now C = min(i/, q;)(2/C^jvs)^^''- 

Second Step. In order to integrate (|3.59p . we need to analyze the exponent of 
nonlinearity fi given by l|3.58p . But, firstly, we recall that /i is written in terms 
of the interpolation exponent 9. According to Lemma Fl. II with p = q = 2 and 
r = a, and also l|3.56p . the admissible range of 9 shows us that < 9 < I iS 
< cr < 2 for AT 1 or iV > 3, and < 6* < iff < cr < 2 for TV = 2. 
In consequence, within these values of cr, we see that ^ > 1 iS 1 < a < 2. The 
other possible values of cr, i.e. < ct < 1 go out of our initial assumption that 
(7 > 1. 

Third Step. An explicit integration of (|3.59p between t = and t lead us to 
E{t) < (e{Q)^ -Ct)~ , C= mm{iy, a) (^7^) " . (3.60) 
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Notice that, from (|2.44p . the initial energy E{0) is finite, and 

n 4+(2-cr)7V 

^'>'^'<ir^i- 2(2-.) ■ ^'-''^ 

Using ^A^, 1338]) and ([3^ . we prove that the right-hand side of 
vanishes for 

_ J 2[(2-<T)W + 2gl 

^ = r ^ "0 m ^t:^ { ^^^^7 ^ ' 3.62 

C 4(2 — cr) •awnyy^a) 

which proves the first assertion with t* given by l|3.62p . 

Fourth Step. If f 7^ 0, we use Schwarz's and Cauchy's inequahties, the later 
with e ~ vl2, to obtain for a.a. t G [0, T] 



\x{t) ■ f{t) dx 



< WumvWmWv < pVu(i)||^(o) + ^llfWIIv'- (s.es) 



Then, from f^A5\i and ^Ml, 

^E{t) + CiE2,At) < C2\\{{t)\\l> fora.a. [0,T], (3.64) 

where Ci = min(i//2,a) and C2 = 1/(2;^). Using p.52p and l|3.57p . we obtain 
the following non-homogeneous ordinary differential inequality 

1 

^E{t) + CiEity^t" < Cse^ (1 - -V ' for a.a. t e [0, T] , (3.65) 
at \ tf J ^ 

where Ci = ™"('^/^'°) g^jj(j (j^ — 1/(2;^). To analyze (I3.65p . we need the 
following result (see the proof in Antontsev et al. [U §1.2.2]). 

Lemma 3.1 Let 5 > such that {tf—d,tf+S) C [0,T] and assume E G 
W^'^ {tf — S,tf+ S) satisfies the differential inequality 



-Eit)+ip{Eit))<F\^l^l--j j a.e. in{tf-5,tf+5), (3.66) 

where <~p is a continuous non- decreasing function such that 

ip(0) = and [ < 00 , (3.67) 

Ja+ ^[-V 

and the function F satisfies, for some k G (0, 1), to 

F{s) < (l ~k)ip{m:{s)) in{Q,tf), (3.68) 

where 

Vk{s)^Ol^{s) and ^k{s) ^ ' 
Then E{t) = for all t > tf. 
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For the relations (|3.65p and l|3.66p to coincide, let us set 

ip{s) = Cis~' and F{s) — C2^^s~ 
Then clearly l|3.67p is satisfied and we have 



Moreover, l|3.68p is satisfied if 



S 



and 77fc(s) = 



e <1 



9i 



{i-k) 



2v min(;^/2, a) 

\ (2-<t)]V + S 



4+ (2 - cr)Ar 



(2-o-)-W + 2ct 
2(2-c7) 



(3.69) 



Second assertion is thus proved and eo is the term defined by the right-hand 
side of ((3:691) . □ 



Remarks 3.1 1. The results established in Theorem \3.1\ still hold for un- 
bounded domains 51 as far as the used Sobolev type inequalities hold (see Lemma UTl]} . 
It is important that the domain is convex and bounded, at least, in one direc- 
tion. 

2. We could also have considered non-homogeneous boundary conditions, say 
ub, on Tt- But then, in order to carry out the results of Theorem \ 3.1\ we 
would have to assume the existence of ats > such that us = for all t >tB 
and E{tB) < co. In the above proof we only would have to replace the time t ~ 
by t = tB and, for the second assertion, to presume that tf> Ib- 

If we assume, in Theorem EH f G l?{Q,T;V{n)), then we must replace 
l(332l) by 



IfWI 



LP(0) 



t \ 



27V 
N + 2' 



N^2, 



where now ^ depends on p, N and a. Indeed, proceeding such as for p.63p and 
using, in addition, (|1.13p with q — p' , p ^ 2 and 6* = 1, we obtain 



u{t) ■ f{t) dx 



<l|uWII 



LP (O) 



\m\\ 



LP(n) 



< 



C\\^^mi.^ia)\\m\knn)<4^^ml^n)+C{£)\m\\lnn)^ 

For N ~ 2, we assume f e L^{0,T;I?{n)) and by using Poincare's inequality 
ifTTsl) . we get 



u(i) • f{t) dx. 



<l|u(t)||L=(n)l|fWllL 



< 
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The results of Theorem 13.11 can be extended to the limit case of cr = 1. In 
fact, if (T = 1 in (|3.58p . then /i = 1 + 2/{N + 2) and all the reasoning of the 
above proof can be applied. On the other hand, if cr = 2 these results are no 
longer vaHd, because, from l|3.62p and p.69p . t* ^ cg and eo ^ as cr ^ 2. In 
this case, taking cr = 2 in p.SSp . we obtain from (|3.59p 

j^E{t) + CE{t) < for a.a. te[0,T], C = 2 mm{iy, a) /C^Ng . (3.70) 

Notice that i?2,o-(0 < E{t) is trivial for cr = 2. Integrating l|3.70p between t = 
and t > 0, we prove the weak solutions of the modified NS problem (|1.7p - (|1.10p . 
with cr = 2 and f = 0, have the following exponential decay 

Mt)\\-LHn) < Cie-^'' for a.a. t>0, (3.71) 

where Ci = ||uo||L2(n) and C2 = 2 min{i^, a) /C^j^g. The following theorem 
shows us that for cr > 2, the weak solutions of the modified NS problem l|1.7p - 
Ijl.lOp . with f = 0, have a power-time decay. 

Theorem 3.2 (Power decay) Assume f— and Uq ^ H, and let u be a weak 
solution of the modified NS problem ii. 7)) -/ O770)) in the sense of Definition \2.1[ 
If a > 2, then there exist positive constants Ci and C2 such that 

Mt)\\LHn) <{Cit + C2r^ for a.a. t>0. 

PROOF. Firstly we observe that using Holder's, Sobolev's (|1.13p and algebraic 
l|1.12p inequalities, we obatin 

er + 2 

\\uml.^,,^ < Cs (^j^ (lu(t)r + |Vu(t)H " for a.a. t e [0,r] , 

(3.72) 

where Cs is the constant resulting from applying Sobolev's inequality (|1.13p . 
Then (|3.53p and l|3.72p lead us to the homogeneous ordinary differential inequal- 
ity 

j^E{t) + CE{t)^ < for a.a. t e [0, T] , (3.73) 

where C = (2/C,s)2'^/('^+2) min(:^, a). Integrating l(3?73l) between t = and i > 
ends the proof, where 

C,= ^—j mm(.,«)— ^ and = ^ j .□ 

Remarks 3.2 1. By virtue of using Sobolev's inequality, ^3. 72]) holds for any 
a > if N > 3 or for any a > 1 if N = 2, and the assumption cr > 2 

satisfies these requirements. 

2. In the limit case a = 2, \3. 7g|) becomes a linear differential inequality anal- 
ogous to fg. yO|) and, again, we obtain the exponential decay fS. 71]) . where we 
only have to replace Cqj^^ by Cs- 
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In Theorem 13. II we have seen that for 1 < a < 2 and f = or for f satisfying 
to (|3.52p it was possible to estabUsh that the weak solutions of the modified 
NS problem l|1.7p - l|1.10p extinct in a finite time. If cr > 2 and f = it was 
possible to prove, in Theorem 13.21 that the weak solutions exponentially decay 
if cr = 2 and decay at a power rate if cr > 2. Therefore it is reasonable to ask 
what happens if f 7^ merely belongs to a suitable function space. In order to 
work this case, let us go back to (|2.45|) . Using there Cauchy's inequality with 
e = 1^/2, we obtain 



dt 



E{t) + CiE2,a{t) < C2 [ \{{t)\^ djc for a.a. t e [0, T] , (3.74) 
Jn 



where E2,a{t) is defined in i|3.54p and now C ~ min(j//2,Q;) and C2 = l/{2u). 
Now we assume that 

/ |f(i)|^ dx < Cf for a.a. t e [0,T] , (3.75) 
Jn 

where Cf is a positive constant. If 1 < cr < 2 we use l|3.57p and if cr > 2 we use 
l|3.72p . In any case, we obtain from l|3.74p and (|3.75p 

j^E{t) + CiE{ty < C2 for a.a. t £ [0, T] , (3.76) 

where 

(2-cr)iV + 2cr , 2cr , , ^ 

7 = 777 TTTT if K < 2 and 7 = if cr > 2 , 3.77) 

' 4+(2-cr)iV ' cr + 2 ^ \ J 

Ci = (^Q^ (2'") ^^'^ ^ ^ cr > 1 , (3.78) 

and where C = C^j^g or C = Cs if 1 < cr < 2 or cr > 2, respectively. Notice 
that in the Hmit case of cr = 2 both expressions of 7 in (|3.77p converge to the 
same value 7 = 2. Let us set now 

^^E{t) + C,E{tr = C2 ^ ^^E{t) = C2 - C,E{tr ■■= A{t) . (3.79) 

If A(f) < or A(t) > at some time t (possibly different), then E{t) is decreasing 
or increasing, respectively, at that time. In consequence, the asymptotically 
stable equilibrium of (|3.79p is reached when 

A(t) = ^ E{t) = (^) ' = C (- .^'^ ^ ] " := ^ = , (3.80) 

^' \Ci) \2vTahi{v/2,a) ) 2 ^ > 

where Ci, C2 and C are given in l|3.78p . A simple analysis of l|3.79p and i|3.80p 
shows us that if there exists a positive time such that Eito) < E^, then 
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< E(t) < E^, for all time t > to and, consequently, E{t) E^ as t ^ oo. In 
this case, we are done and we obtain for all t > 

4 + (2~g)W 



or 



0- + 2 

\\nmhin)<^Cs(- r^^Sr^V" if cT>2. 

The reciprocal case shall be studied in the next theorem. 

Theorem 3.3 (Exponential decay) Assume uq ^ H and f ^ satisfies to 
113.75]) . Let u be a weak solution of the modified NS problem ( ti. 7)) - f7T75)) in the 
sense of Definition \2.1l In addition assume that there exists a positive time to 



such that ||u(to)||^2(fj-) > £*■ Then there exists a positive constant C such that 

||«(*)ll^(o) < (ll«(io)|li2(o) -f,)e-^(*-*«'+£, for all t> to, (3.81) 
where is given in ^3.80]) . 

PROOF. First Step. Again a simple analysis of (|3.79p and (|3.80p shows us that 
if there exists a positive time to such that E{to) > then E{t) > E^ for all 
time t > to and, consequently, E{t) \ E'* as t oo. In order to simplify the 
notations, let us set E{t) = E, -f^Eit) = E' and f{E) CiE'^. First, we 
observe that with these notations we have from l|3.79p and p.80p . 

E' + CiE'' = C2 ^ {E-E,y + Ci— -^{E-E,)=0. (3.82) 

E — E^, 

Second Step. Using (|3.80p . we can prove that 

C,^-^^l^^l-^>nE.)^^C^C^ iff 7>1. (3.83) 

According to the expressions of 7 (see l|3.77p ). we can see that p.83p holds iff 
a > 2. Then from (|3.82p and l|3.83p we derive the following linear differential 
inequality 

{E - E^y + C{E -E^) <0, C = 7 C^C^ . (3.84) 

Integrating (j3.84p between to and t > to, we obtain l|3.8ip where, from l|3.78p 
and ([3^ . 

2a 2 . /i^ \^ { Cf\ ^ 



2Cs V2' / \2v 

Third Step. For 1 < a < 2 and in addition to what we have done in the First 
Step, let us set i?(to) = Eo- We then prove 

^^i;^>^^lh^ iff 0<7<1, (3.86) 
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and according to l|3.77p . we see that this holds iff 1 < cr < 2. Then from l|3.82p 
and (|3.86p . we derive the following linear differential inequality 

{E-E,y + C{E-E,)<0, C = Ci ^J~ ^* . (3.87) 

Eq — Elf 

Integrating l|3.87p between to and t > to leads us to p. Sip , where from l|3.78p 
and p.87p . the constant C is given by 

(2-<tUV + 2<t (2-g)W + 2g (2-<r)W + 2g 

4+(2-cr)JV 771 4 + (2-<t)N 



^GNsJ V2' / Eq — E. 



C= 7^ min .□ (3.88) 



To conclude this section, we analyze the results obtained in this section in 
terms of a and ly. Theorem 13.11 is the most restrictive, because if one and only 
one of the constants a and v tend to zero, we obtain, from p.62p and p.69p . 

^ oo and e — > and those results fail. But if only a tend to zero, the results 
corresponding to exponential and power-time decays (Theorems 13.21 and 13. 3p 
still hold. This situation corresponds to consider the classical Navier-Stokes 
problem. If tends to zero, the situation gets worse, because we can no longer 
use Sobolev-type inequalities. 



4 Conclusions 

Throughout this paper we have seen that the introduction of the absorption 
term |u|°'~^u in the momentum equation of the Navier-Stkes problem allows us 
to obtain different time properties for the weak solutions. If 1 < cr < 2 we have 
improved the known results for the NS problem and we obtained an extinction 
in a finite time, whereas for cr > 2 we have obtained exponential or power-time 
decays. The problem lies in how to justify, from the physical point of view, the 
appearance of this term in the momentum equation. We may think of it as ex- 
pressing somehow a sink that dissipates kinetics energy. Or it can correspond to 
a physical body inside the fiuid, or at the boundary, which makes the fiow slower 
when the time is passing by or even stop it in a finite time. Other possibility 
is to consider the absorption term as resulting from a real or fictitious force as 
already mentioned in (jl.lip . For instance, in Geophysical flows the flctitious 
CorioHs force resulting from the Coriolis acceleration is given by 2r2xu, where 
CI is the Earth angular velocity, and a|u|'^~^u approximates the CorioHs force 
if cr = 2. In this case, we only obtain exponential and power-time decays. From 
the theoretical point of view, the results of this paper can be extended to a great 
variety of Fluid Mechanics problems modifled by introducing in the momentum 
equation an absorption like term. In OHveira [20] we consider the Oberbeck- 
Boussinesq problem modifled by the thermo-absorption term q;|u|'^(^^~^u, where 
9 is the absolute temperature. Here cr is a temperature depending function with 
Lipschitz regularity and such that a{0) > 1 for all 6* S R. Under the assumption 
that 1 < cr^ < a{9) < cr+ < 2 for all 6* G M, where a" and cr+ are constants, 
we are able to establish the same asymptotic stability properties. We may also 
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consider the non-homogeneous Navier-Stokes problem modified by an absorp- 
tion term. In this case, we only need to assume that the density function p 
is bounded: ^ < p , Pn < Cp, where Cp is a positive constant and po is the 
initial density. These properties extend also for the modified NS problem sup- 
plemented with the slip boundary conditions: u • n = and u • r = /3~^t • r 
on Ft, where n and t denote, respectively, unit normal and tangential vectors 
to the boundary Sfi, t = n • T is the stress vector and /? is a coefficient with no 
defined sign (see Antontsev and Oliveira [2]). As far as the behavior in space 
of the weak solutions of such modified problems, at the moment, we are able 
to prove analogous properties only for the 2D stationary modified NS problem 
(see Antontsev and Oliveira [3]). 
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